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Abstract 

The  self-discharge  rate  of  aerospace  flight-quality  nickel-hydrogen  batteries  is  limited  by  hydrogen  diffusion  within  the  nickel  electrode 
active  material.  A  diffusion-limited  reaction  model  is  developed  which  accounts  for  the  observed  self-discharge  behavior.  Effective  hydrogen 
diffusion  coefficients  were  calculated  from  the  open-circuit  self-discharge  data  on  Hubble  space  telescope  flight-qualified  nickel-hydrogen 
cells  and  ranged  from  1 . 1 2  X 10' 10  cm2/sec  at  25  °C  to  7.45  X 10' 12  cm2/sec  at  0  °C. 
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1.  Introduction 

Historically,  the  self-discharge  reaction  rate  of  nickel- 
hydrogen  batteries  has  been  proposed  to  be  under  kinetic  or 
electrochemical  control  [  1-8] .  However,  calorimetric  meas¬ 
urements  on  charged  nickel  electrodes  under  constant  hydro¬ 
gen  pressure  indicates  that  the  self-discharge  reaction  rate  is 
not  constant  with  time  [6,9-11].  If  the  first-order  kinetics 
with  respect  to  hydrogen  were  controlling  the  self-discharge, 
the  reaction  rate  would  be  constant  with  time  under  constant 
hydrogen  pressure.  Second-order  overall  kinetics  [4,6],  first- 
order  kinetics  with  respect  to  hydrogen  and  p-NiOOH,  does 
not  account  for  the  effect  of  hydrogen  diffusion  in  the  solid 
nickel  active  material.  In  a  previous  electrochemical  model 
[7,8],  the  effect  of  non-catalytic  chemical  reaction  and  dif¬ 
fusion  has  also  not  been  accounted  for.  A  simple  chemical 
reaction  and  diffusion  model  to  explain  the  self-discharge 
phenomenon  has  been  previously  published  [12].  An 
updated  mathematical  description  of  hydrogen  diffusion  and 
non-catalytic  chemical  reaction  within  the  nickel  electrode 
active  material  to  describe  the  self-dischargeof  nickel-hydro¬ 
gen  constant  volume  batteries  on  open-circuit  at  various  tem¬ 
peratures  is  presented  in  this  paper.  Self-discharge  data  and 
their  correlation  with  the  simple  diffusion-limited  reaction 
model  are  also  presented.  The  self-discharge  data  are  from 
the  open  literature  on  a  constant  hydrogen  pressure  nickel- 
hydrogen  cell  and  on  nickel-hydrogen  batteries  that  are  in 
use  on  the  Hubble  Space  Telescope  [5,6,10] . 


2.  Self-disc  harge  model 

To  promote  rapid  oxygen  recombination  on  overcharge, 
the  nickel  electrode  of  the  nickel-hydrogen  battery  operates 
in  an  ‘electrolyte-starved’  condition  [  13,14].  Therefore,  the 
electrolyte  is  only  present  in  the  micropores  of  the  active 
material  and  the  macropores  of  the  nickel  electrode  are  free 
of  electrolyte  for  rapid  transport  of  oxygen  produced  on  over¬ 
charge.  Assuming  that  the  nickel  active  material /electrolyte 
system  is  a  homogeneous  mixture,  the  nickel  electrode  can 
be  modeled  as  shown  in  Fig.  1.  Due  to  the  low  mass  transfer 
rate,  pseudo  steady-state  thermodynamic  equilibrium  is 
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Fig.  1 .  Schematic  diagram  of  nickel  electrode. 
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assumed  for  hydrogen  and  the  concentration  of  the  dissolved 
gas  in  the  electrolyte  at  the  interface  (jr  =  0)  is  represented 
by  Henry’s  law.  Hydrogen  is  dissolved  in  the  electrolyte  and 
transported  in  the  homogeneous  active  material /electrolyte 
mixture  and  chemically  reacts  with  fi-NiOOH  to  g-Ni  ( OH )  2. 
This  reaction  is  given  by  Eq.  ( 1 ): 

H2  +  2(3-NiOOH  ->  2($-Ni  (OH  )2  (1) 


A  one-dimensional  mass  balance  on  hydrogen  at  x= *,  (Fig. 
1),  with  constant  density,  porosity,  and  diffusivity  and  with 
negligible  convective  transport,  results  in: 


ac. 


(2) 


A  mass  balance  on  the  nickel  oxyhydroxide  yields: 

3C, 

-e=-*2=-  2K,  (3) 

ot 


In  dimensionless  form,  Eqs.  (2)  and  (3)  become: 

(l£iA9A=ltii  d* 

\  C10  )  dr  M2  9f2  K' 


(4) 


Cl  =  C,s(‘"3 

Eq.  (8)  reduces  to: 

I  d«,^  SDcCls 
l  d  t  j  x, 

The  consumption  rale  of  the  nickel  oxyhydroxide  yields: 

(-2K-® 

Substituting  Eq.  ( 10)  into  Eq.  ( 1 1 )  yields: 

— =2^^  (12) 

Qt  p2X{ 

The  concentration  of  hydrogen  in  the  electrolyte  at  the  inter¬ 
face  (jt=0)  is  given  by  Henry’s  law  as  shown  below: 

P^HCU  (13) 

From  geometrical  considerations,  the  fractional  conversion 
of  3-NiOOH  is  given  by: 


(9) 

(10) 


Wi 

3t 


--i?r 


(5) 


When  1,  the  hydrogen  diffusion  rate  is  much  faster 
than  the  kinetic  reaction  rate  and  the  homogeneous  model 
may  be  applied.  When  A/2»  1,  the  kinetic  reaction  rate  is 
much  faster  than  the  diffusion  rate  and  the  overall  reaction 
rate  is  therefore  limited  by  diffusion.  When  the  diffusion  and 
kinetic  reaction  rates  are  of  the  same  order  of  magnitude 
(Af2~l),  then  Eqs.  (4)  and  (5)  must  be  simultaneously 
solved  in  their  full  form  [  15] . 

From  calorimetric  measurements  on  ‘starved-electrolyte’ 
nickel  electrodes,  M2  appears  to  be  much  greater  than  one 
[6,9].  Thus,  the  pseudo  steady-state  approximation  of  the 
hydrogen  mass  balance  leads  to: 


tii  jpc, 

df2  dr 


(6) 


Integration  of  Eq.  (6)  results  in  the  pseudo  steady-state 
approximation  of  a  constant  flux  with  respect  to  the  spatial 
jc-coordinate. 


N  i  = 


= constant 


(7) 


The  instantaneous  consumption  rate  of  hydrogen  is  then  given 
by: 


(14) 


Substituting  for  C,  s  and  x{  from  Eqs.  ( 1 3)  and  ( 14),  respec¬ 
tively,  into  Eq.  ( 12)  yields: 


dX2  jg/l 
d  /  ~Hp2X2Ll 


(15) 


For  constant  hydrogen  pressure  systems,  integration  of 
Eq.  (15)  leads  to  the  unconverted  fraction  of  nickel  oxy¬ 
hydroxide  as  shown  below: 


(16) 


For  a  constant  hydrogen  pressure  battery  of  an  initial  fully 
charged  capacity  of  A0,  the  capacity  with  time,  A,  is  then 
given  by: 

<>’> 


Thus,  a  plot  of  remaining  capacity  versus  the  square  root  of 
open-circuit  self-discharge  time  would  yield  a  straight  line. 

For  constant  volume  batteries,  the  fractional  conversion  of 
hydrogen  and  the  nickel  oxyhydroxide  may  be  defined  as: 


X. 


Thus: 


(IB) 


Integration  of  Eq.  (6)  from  a=0  to  jc=jc,  and  noting  that 
C |  =C,  „  at  .r=0,  and  Cj  =  0  at  respectively,  yields 
the  following  equation  for  the  hydrogen  concentration 
distribution. 


Substituting  Eqs.  ( 18)  and  ( 19)  into  Eq.  ( 15)  yields: 
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fJYo -P,V  dP,\  2PC(PI„-P,I,)- 

l  P,  A  dfj 


(20) 


Integration  of  Eq.  (20)  from/ =0  tot  and  noting  that  P,  =  PU) 
at  f = 0  and  Pj  =  P,  at  t  - 1,  yields  the  following  equation: 


\Pj  \  P,. 0  I  \ 


(21) 


Thus,  a  plot  of  the  left-hand  side  of  Eq.  (21)  versus  time 
should  yield  a  straight  line. 

Eq.  ( 2 1 )  can  be  further  simplified  by  performing  a  Taylor 
expansion  of  the  natural  log  portion  of  the  left-hand  side  and 
retaining  the  first  three  terms.  Thus: 


Substituting  Eq.  (22)  intoEq.  (21)  and  simplifying  yields: 


P,=P,.  o- 


l4Pc(Pl.o-Pi.p)?J>,, 

HfhLl 


{> 


(23) 


Thus,  for  open-circuit  pressures  not  too  far  from  full  charge, 
the  self  discharge  behavior  of  constant-volume  batteries  can 
be  approximated  by  Eq.  (23).  Note  that  the  only  variable 
term  on  the  right-hand  side  of  Eq.  (23)  is  time. 

Using  Eqs,  (18)  and  (23),  the  unconverted  fraction  of 
p-NiOOH,  here  defined  as  the  capacity  retained  (CP),  is 
given  by: 

=  (24) 

V  HpA 


Thus,  Eq.  (24)  can  be  used  to  approximate  the  self-discharg¬ 
ing  behavior  of  a  constant- volume  nickel-hydrogen  battery 
on  open-circuit  starting  from  a  full  charge. 


3.  Validation  of  the  diffusion-limited  reaction  mode! 

To  verify  the  proposed  mathematical  model,  experimental 
data  from  the  open  literature  on  self-discharge  of  nickel- 
hydrogen  batteries  were  obtained.  Data  on  constant  pressure 
and  volume  systems  were  used  to  generate  Fig.  2  and  Fig.  3, 
respectively  [5,10], 

The  set  of  data  for  Fig.  2  was  obtained  from  a  nickel- 
hydrogen  cell  with  31  wt  %  KOH  in  water  as  the  electrolyte, 
at  a  constant  hydrogen  pressure  of  1 50  psig  and  a  temperature 
of  25  °C.  The  diffusion-limited  reaction  model  behavior  for 
a  constant  pressure  battery  is  given  by  Eq.  (17).  Using  a 
linear  least-squares  regression,  the  model  equation  was  fitted 
with  a  coefficient  of  determination  of  0.998. 

Data  sets  for  Fig.  3  were  obtained  from  the  nickel-hydro¬ 
gen  cells  that  passed  flight  qualification  testing  for  the  Hubble 
space  telescope.  The  cells  used  in  these  tests  were  from  the 
same  cell  lots  as  the  cells  that  are  presently  in  orbit.  To  obtain 
the  0  and  25  °C  open-circuit  data,  the  cells  were  charged  at  0 
°C.  The  20  °C  open-circuit  data  were  obtained  by  charging 
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Fig.  2.  Nickel-hydrogen  cell  self-discharge  at  constant  pressure  of  150  psig. 


Fig.  X  Hubble  space  telescope  nickel-hydrogen  batteries  charge  retention 
test. 

the  cells  at  10  °C.  Approximation  to  the  diffusion-limited 
reaction  model  behavior  for  a  constant  volume  battery  is 
given  by  Eq.  (24).  Using  linear  least-squares  regression,  the 
model  equation  was  fitted  to  the  data  with  coefficients  of 
determination  of  0.994, 0.997,  and  0.999  for  0, 20,  and  25  °C 
open-circuit  data,  respectively. 

For  a  Henry’s  law  constant  of  8  9 1 8  725  atm  cm3/ (g  mol) 
[  16],  nickel  oxyhydroxide  molar  concentration  of  0.042  g 
mol/cm3,  active  material  thickness  of  5X  IQ"4  cm,  and  a 
pressure  of  1 1 .2  atm,  the  calculated  effective  diffusivity  from 
the  slope  of  Fig.  2  is  2.3  X 10" 10  cm2/sec.  For  a  full  charge 
pressure  of  approximately  68  atm  the  calculated  effective 
diffusivities  from  the  slopes  of  Fig.  3  are  1.12X10"10, 
9.67  X 10- 1 and  7.45  X 1 0" 12  cm2/sec  for  25, 20,  and  0  °C, 
respectively. 


4.  Discussion  and  conclusions 


The  magnitude  of  the  hydrogen  diffusivity  calculated  from 
the  slopes  of  Fig.  2  and  Fig.  3  is  similar  to  other  investigators 
values  who  used  different  techniques  [  17-19].  Thus,  a  sim¬ 
ple  diffusion-limited  reaction  model  adequately  explains  the 
self-discharge  phenomena  in  nickel-hydrogen  batteries. 

The  use  of  cobalt  hydroxide  in  the  nickel  electrodes  for 
these  flight-qualified  batteries  significantly  reduces  the 
amount  of  oxygen  generated  near  the  end  of  charge,  thus 
neglecting  the  eff  of  oxygen  in  the  model  is  justified. 
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Since  the  self-discharge  rate  of  nickel-hydrogen  batteries 
is  controlled  by  the  diffusion  rate  of  hydrogen  in  the  electrode 
material,  there  are  several  ways  to  minimize  this  self-dis¬ 
charge  rate.  They  are  ( i )  increasing  the  active  material  thick¬ 
ness  on  the  nickel  sinter  (ii)  decreasing  hydrogen  operating 
pressure:  (iii)  decreasing  battery  operating  temperature,  and 
(iv)  using  a  hydrogen-limited  cell  design. 


5.  List  of  symbols 

A  ampere  hour  capacity 

A{)  initial  ampere  hour  capacity 

C,  hydrogen  concentration 

C2  p-NiOOH  concentration 

Cl>0  initial  hydrogen  gas  concentration 

CJtP  precharge  hydrogen  gas  concentration 

C,.s  hydrogen  concentration  in  the  electrolyte  at  Jt  =  0 
(see  Fig.  1) 

^2,0  initial  P-NiOOH  concentration 

CR  capacity  retained  ( fraction  unconverted ) 

Dc  effective  diffusivity 

H  Henry’s  law  constant 

first  order  rate  constant 
L.a  thickness  of  active  material 

M2  #i(Ci,s;  C2{))I(DcC]  jLl)\  the  Thiele  modulus 

squared  whi^h  is  indicative  of  the  ratio  of  the 
maximum  reaction  rate  over  a  characteristic 
diffusion  rate 

n,  moles  of  hydrogen 

jV,  hydrogen  flux  rate 

«2  moles  of  p-NiOOH 

P !  absolute  hydrogen  pressure  corrected  for 

compressibility 

Pl  Q  initial  corrected  hydrogen  pressure 

P,  p  precharge  corrected  hydrogen  pressure 

/?,  rate  of  hydrogen  consumption  per  unit  total 
volume 

Ri  *  RJRiiCx^  C2t0);  dimensionless  hydrogen 
consumption  rate 

R2  rate  of  P-NiOOH  consumption 

S  total  ( grai  n  +  pore )  planar  surface  area  of  the 
nickel  active  material  perpendicular  to  the  jc- 
coordinate 
t  time 

x  distance  in  nickel  active  material 


x{  distance  in  nickel  active  material  where  diffusion- 
limited  reaction  front  occurs 
X2  fractional  conversion  of  p-NiOOH 

Vb  battery  gas  volume 

Greek  letters 

e  p-Ni(  OH)  2  porosity 

</>,  C,  /C|  S;  dimensionless  hydrogen  concentration 

cf) 2  C:/C2,)i  dimensionless  p-NiOOH  concentration 

£  x/La\  dimensionless  distance 

p2  p-NiOOH  bulk  molar  density 

t  2/?j(Cj  s;  C2,))//C2n;  dimensionless  time 
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